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THE  MECHANICS  OF  THE  TRIAXIAL  TEST  FOR  SOIL S* 


R.  M.  HAYTHORNTHWAITE** 


Abstract 

Stress  distributions  at  failure  are  not  uniquely  defined  in  terms  of 
equilibrium  and  a  failure  criterion  such  as  the  Coulomb  law  which  is  commonly 
used  for  soils.  This  fact  is  illustrated  in  the  paper  by  constructing  alterna¬ 
tive,  non-uniform  stress  distributions  for  the  triaxial  test  which  lead  to 
different  apparent  angles  of  friction.  A  suitable  flow  rule  is  introduced  and 
complete  solutions  are  obtained  in  which  the  external  pressures  causing  failure 
are  uniquely  defined.  Several  alternative  velocity  fields  are  developed  and  the 
extent  of  the  deformable  zone  is  established.  The  analysis  is  used  to  discuss 
the  significance  of  tests  on  sand  where  the  specimens  are  caused  to  fail  with 
axial  extension  rather  than  with  the  more  conventional  axial  compression. 


This  paper  reports  work  sponsored  by  the  Land  Locomotion  Research  Branch, 
Ordnance  Tank-Automotive  Command,  Detroit  Arsenal,  under  Contract 
DA-19-020-0RD-U566 . 


Associate  Professor,  Division  of  Engineering,  Brown  University, 
Providence  12,  R.  I. 


Introduction 


It  is  customary  to  base  both  the  analysis  of  the  stability  of  soil 
masses  and  the  interpretation  of  soil  tests  upon  the  Coulomb  theory  of  internal 
friction^  This  theory  provides  what  in  the  terminology  of  plasticity  theory  is 
called  a  yield  or  failure  surface  and,  together  with  the  equilibrium  conditions, 
it  enables  the  failure  loads  to  be  confuted  whenever  the  problem  is  statically 
determinate.  Other  problems  which  are  not  intrinsically  so  can  be  rendered 
statically  determinate  by  invoking  the  hypothesis  of  Haar  and  v.Karman^,  which 
asserts  that  the  intermediate  principal  stress  is  equal  to  either  the  largest  or 
the  smallest  principal  stress.  This  hypothesis  is  particularly  useful  in  axially 
symmetric  problems. 

The  conventional  approach  outlined  above  suffers  from  the  drawback  that 
the  solutions  so  obtained  are  not  necessarily  unique.  To  ensure  uniqueness,  a 
stress-strain  relation  of  a  particular  type  is  required  and  the  establishment  of 
such  a  relation,  if  it  exists,  is  one  of  the  most  important  of  the  unsolved  pro¬ 
blems  of  soil  mechanics. 

This  paper  is  concerned  with  an  examination  of  the  so-called  triaxial 
test  in  the  light  of  a  particular  stress-strain  relation  which  would  ensure  both 
stability  and  uniqueness  of  the  stresses  if  its  validity  could  once  be  established 
by  means  of  tests.  The  analysis  is  used  to  discuss  the  significance  of  tests  on 
sand  in  which  the  specimens  are  caused  to  fail  with  axial  extension  rather  than 
with  the  more  usual  axial  compression. 

The  triaxial  test.  A  detailed  account  of  the  triaxial  test  has  been  presented  by 
Bishop  and  Henkel ^  in  a  recent  monograph.  The  essentials  of  the  device  are 
indicated  in  Fig.  1.  A  cylindrical  specimen  of  soil  is  enclosed  within  a  rubber 
membrane  and  mounted  between  stiff  end  plattens ,  Hydrostatic  pressure  is  applied 
by  means  of  an  enveloping  fluid,  which  is  isolated  from  the  specimen  by  a  rubber 
membrane,  and  an  extra  axial  thrust  or  pull  can  be  exerted  on  the  end  plattens. 
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During  a  test,  the  relative  displacement  of  the  end  plattens  is  observed,  and 
sometimes  the  change  in  volume  of  the  specimen  and  its  change  in  diameter  at  a 
particular  height  are  also  observed.  Testing  proceeds  by  varying  the  axial  load 
or  the  lateral  pressure  or,  quite  commonly,  by  imposing  axial  displacement  of  a 
platten  and  observing  the  resultant  axial  loads. 

This  apparatus  now  occupies  a  central  position  among  the  devices  used 
to  measure  the  strength  of  soils.  It  appears  to  have  developed  from  a  machine 
designed  at  the  Prussian  Waterways  Experimental  Station  for  the  purpose  of  study¬ 
ing  the  consolidation  of  clays  under  conditions  of  negligible  side  friction^ -L 
In  this  first  apparatus  the  surrounding  liquid  was  entirely  confined,  and  tempera¬ 
tures  and  leakage  had  to  be  closely  controlled  to  obtain  consistent  results. 
Several  investigators  recognized  that  the  apparatus  could  be  used  to  measure  the 
ratio  of  the  axial  and  lateral  pressures  both  prior  to  and  at  failure,  the  first 
results  being  apparently  those  published  by  Stanton  and  Hveem^]  in  193J4.  Posi¬ 
tive  control  over  the  lateral  pressure  was  developed  independently  by  Rendulic  ^ 

£6] 

working  in  Berlin,  and  by  Housel  at  the  University  of  Michigan, 

Throughout  the  literature,  it  is  assumed  that  a  homogeneous  state  of 
stress  is  produced  in  the  triaxial  soil  specimen.  This  point  of  view  appears  to 
have  been  taken  over  from  the  earlier  studies  of  rocks  by  von  Karman  and  others 
and  does  not  appear  to  have  been  questioned  since  that  time.  Strictly  speaking, 
however,  the  problem  posed  by  the  test  situation  is  statically  indeterminate  be¬ 
cause  only  the  total  thrust  on  the  end  plattens,  not  the  pressure  distribution, 
is  given  as  a  boundary  condition.  Even  allowing  circular  symmetry  and  similarity 
of  conditions  on  every  normal  cross-section,  there  are  still  three  unknown  stresses 
or  ,  Og  and  oz  (see  Fig,  2)  and  only  two  equilibrium  equations: 


dar  +  CTr-  ~  g9 
dr  r 

^z 

dz 


( 

0  I 

\ 


(1) 


3 


If  now  oq  is  equated  arbitrarily  to  or  ,  or  is  then  constant  by  the  first  of 
(l)  and  the  stress  is  homogeneous.  This  step  has  never  been  established  a  priori , 
but  it  has  seemed  reasonable  in  view  of  its  correctness  in  the  elastic  range 
(at  sections  remote  from  the  ends  of  the  specimen).  Moreover,  some  such  assump¬ 
tion  has  been  necessary  before  information  about  yield  could  be  extracted  from 
the  test. 

The  Coulomb  yield  criterion.  The  weight  of  evidence  from  the  triaxial  test 

[71 

supports  the  friction  criterion  of  Coulomb  ,  although  minor  departures  are 
common.  In  this  criterion  the  shear  stress  1 0  causing  slip  on  any  plane  is  taken 
as  the  sum  of  a  constant  value,  termed  the  cohesion,  and  an  additional  amount 
which  is  proportional  to  the  normal  pressure  acting  across  the  plane;  thus 

t:0  =  c  -  a  tan  <p  (2) 

where  c  is  the  cohesion  and  o  the  tensile  stress  across  the  plane.  The  angle  <p 

is  called  the  angle  of  friction  of  the  material.  This  relationship  follows  a 

straight  line  in  the  a  -  t  plane  (Fig.  3),  in  which  the  representation  of  stress 

[8] 

at  a  point  due  to  Mohr  can  also  be  drawn.  A  stress  state  for  which  yield  is 
incipient  on  some  plane  will  be  represented  by  a  circle  touching  the  failure  line, 
such  as  that  shown  in  the  Figure .  From  the  geometry  of  the  triangle  ABC , 

(a3”a1  )/2  *  [c  cot<P-(a1+a3)/2]  sin  cp  where  o ^  ,  a3  are  principal  stresses; 
hence  after  rearrangement 

Oy  «  Cretan2  (^+|)  -  2c  tan  (  ^  +-|.)  .  (3) 

It  is  a  property  of  the  Mohr  circle  that  the  angle  subtended  at  B,  Fig. 3, 
represents  twice  the  angle  of  inclination  of  the  plane  on  which  the  corresponding 
stresses  act.  Thus  the  critical  condition  represented  by  point  C  in  Fig,  3  occurs 
on  a  plane  the  normal  of  which  subtends  an  angle  (  ^  )  with  the  normal  to  the 

plane  on  which  <7^  acts. 
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When  the  stresses  in  the  triaxial  specimen  are  assumed  to  be  homogeneous , 
the  lateral  and  axial  pressures  can  be  related  at  once  through  (3).  It  will  be 
shown  in  the  section  on  Limiting  Equilibrium  States  that  this  solution  is  not 
unique . 

Stable  plastic  flow.  Lack  of  uniqueness  can  be  resolved  by  introducing 
a  stress-strain  relation.  It  is  not  known  whether  a  stress-strain  relation  exists 
for  any  particular  soil:  the  purpose  here  is  to  make  an  initial  assumption  and 
then  to  examine  the  consequences,  with  a  view  to  obtaining  predictions  that  can 
be  compared  with  experiments. 

We  shall  assume  an  isotropic ,  stable  soil  which  yields  at  stress  com¬ 
binations  that  are  constant  and  not  a  function  of  the  strains  or  of  the  stress 
history.  This  is  the  specification  for  a  so-called  ! ideally  plastic®  material. 

By  a  stable  material  is  meant  one  from  which  work  cannot  be  extracted  in  any 
loading  program,  and  it  follows  from  this  definition1"^  that  the  yield  surface 
of  the  material  is  convex  when  drawn  in  a  space  with  the  stress  components  as 
the  axes.  Also  any  plastic  strain  increment  vector  superimposed  on  the  stress 
space  with  corresponding  axes  matching  will  lie  in  the  direction  of  the  outwards- 
drawn  normal  to  the  yield  surface.  If  the  yield  surface  has  corners,  the  vector 
may  lie  between  adjacent  normals. 

When  the  yield  surface  has  a  continuously  turning  tangent,  the  stresses 

in  a  deforming  body  are  uniquely  defined  .  If  it  contains  flats,  there  may 

be  some  latitude,  although  the  choice  may  often  be  restricted  by  the  requirements 
[191 

of  equilibrium  .  For  a  body  with  specified  surface  velocities  which  are  uni¬ 
form,  as  in  the  case  of  the  triaxial  specimen  deforming  between  rigid  plattens, 
it  may  be  shown  that  the  external  tractions  are  unique  despite  the  ambiguity  in 
the  internal  stressing  (see  the  Appendix). 

The  nature  of  the  ideally  plastic  material  which  follows  the  Coulomb 

[10] 

yield  criterion  has  been  investigated  by  Shield  and  certain  results  are 
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quoted  here  for  reference. 

Figure  U  shows  the  intersection  of  the  yield  surface  (3)  with  a  plane 

cr^  =  const,,  where  is  one  of  the  principal  stresses  (not  necessarily  the 

intermediate  principal  stress).  The  regions  corresponding  to  various  orderings 

of  a  ,  <j  and  a  are  indicated.  The  strain  rate  vectors  will  lie  in  the 
a  *  b  c 

directions  of  the  outward  drawn  normals  to  the  sides  of  the  cross-section 


shown  provided  the  components  are  related  in  the  manner  indicated  in  Table  I.  In 
this  Table,  the  values  X^,  X^ 

In  each  line  of  Table  I,  the  sum  of  the  strain  rates  is  proportional  to 
N2  -  1  and  the  sum  of  the  absolute  values  to  N2  +  1 j  hence 


are  constants  and  N  =  tan 


6a  +  %  +  ec  "  (I  6  J  +  I  eb  I  *  *  ec  I  )  sin  <P  0  ♦ 

At  least  one  of  the  three  components  is  positive  and  at  least  one  negative. 
Denote  the  positive  component  by  ea  and  the  negative  by  e^.  When  the  third 
coirqponent  is  positive 

6a  +  ep  +  eY  tan2  (j~  ♦  =  0  (U) 

and  when  is  negative, 

ea  tan2  ^  f  _  |  ^  +  ep  +  6r  "  0  • 

At  the  apex,  the  rate  of  energy  dissipation  D  *  is 


D  **  c  cot  <p  ( e-j  +  62  +  63  )  ■  c  6  cot  <p  (6) 

where  6  is  the  dilatation  rate  per  unit  volume.  Expression  (6)  also  applies  at 
all  other  points  on  the  pyramid  because  the  projection  of  the  stresses  in  direc¬ 
tions  normal  to  the  sides  is  always  the  same. 

Using  the  ideally  plastic  material  described  above,  it  is  possible  to 
obtain  a  unique  solution  for  the  applied  pressures  in  the  triaxial  test.  This 
will  be  done  below,  but  first  the  lack  of  uniqueness  of  the  so-called  'limiting 
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equilibrium  states*  will  be  demonstrated  by  exhibiting  alternative  solutions „ 


Limiting  Equilibrium  States 


A  limiting  equilibrium  state  will  be  defined  as  an  equilibrium  stress 
distribution  in  which  the  soil  reaches  failure  at  a  sufficient  number  of  points 
for  slip  to  appear  possible.  In  the  case  of  the  triaxial  test5  limiting  equili¬ 
brium  states  can  be  based  on  a  homogeneous  stress  distribution*  but  there  are 
alternatives *  based  on  non-homogeneous  distributions *  which  meet  the  hypothesis 
of  Haar  and  v,  Karman  and  a  still  wider  class  of  solutions  that  does  not. 

For  the  purposes  of  illustration  a  general  class  of  solutions  is  devel¬ 
oped  below  in  which  the  principal  stresses  remain  in  constant  ratio  at  all  points 
in  the  specimen*  although  the  stress  level  may  vary  from  point  to  point  in  the 
radial  direction. 


Failure  with  axial 


Shortening  of  the  specimen  will  occur 


when  the  average  axial  pressure  q  =  -cl  exceeds  the  lateral  pressure  p  =  -[cl,]  < 

r  r  =a 

Assuming  >  o  *  (3)  becomes 


ar  =  °z  N  “  2cN 


where  N  =  tan  ( ^  ♦  & ) „  The  position  of  the  intermediate  principal  stress  Oq 
can  be  expressed  by  the  coefficient  k  where 

c9  *  k  ar  +  (l-k)  %  (8) 

and  1  5*  k  >0»  Substituting  (7)  and  (8)  in  the  first  of  (l) *  integrating  and 
making  use  of  the  boundary  condition  [a  ]  _  =  -p  s 

A  I®  s  SI 


a  *  c  cot  <p  -  (p  +  c  cot  (p  )(«£  ) 
x  a 


(N2-l)(k-l) 


oz  =  c  cot  <p  -  N2(p  +  c  cot<p)(^) 


(W2-l) (k-l) 
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It  is  confirmed  that  ar  >  °z  s  because  N2>1.  On  integrating  az  over  the  normal 


cross-section,  the  expression  for  the  mean  axial  pressure  q  becomes 

_  _  2  m2p  ■>  (i*k)sc]  (u) 

(l-k)N2  +  1  *  k 

Setting  k“l(cT0  =  ar),  the  stress  distribution  is  homogeneous  and  (8)  reduces 
to  the  well-known  result 

q  ■  P  tan2(^  +  -|-)  +  2c  tan(J  +  |)  (12) 

while  on  setting  k  ■  0  (aQ  »  a  ), 

q  ■  (l  +  sin  <p)  p  +■  c  cos  cp  (13) 

The  stress  states  giving  rise  to  (12)  and  (13)  both  meet  the  hypothesis  of  Haar 
and  v.  Karman,  in  that  the  intermediate  principal  stress  is  in  each  case  equal 
to  one  of  the  extreme  stresses.  In  addition,  any  number  of  alternative  expres¬ 
sions  may  be  obtained  from  the  stress  distributions  associated  with  other  values 


of  k. 


Failure  with  axial  extension.  Analogous  results  can  be  obtained  for 


the  case  where  the  lateral  pressure  p  exceeds  the  average  axial  pressure  and 
axial  extension  occurs  during  failure.  In  this  case  it  is  assumed  that  crr  <  oz 


and  (3)  becomes 


n  -  2  cN 


Integrating  the  first  of  (1)  after  substituting  (lU), 


c  cot 


(  *  +  wr/jjZ 

<p  -  (p  +  c  cot  <P  )(£)  w 


.  1  /  wrx(^  -  l)(l-k) 

°z  “  c  cot(p-  p  (p  +•  c  cot  <p)(|)  N2 


Here  pj  <  1,  so  otp  <  az.  On  integration  of  (16)  over  the  cross-section 
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(1?) 


_  _  2  [p  -  (l+k)Nc] 
q  (l+k)N2  +  1-  k 

When  k  *  1(oq  -  °r) »  (17)  gives  the  usual  formula  for  the  homogeneous  stress 
distributions  s 

q  =  p  tan2(^  -  |)  -  2c  tan(jj  -  |)  (18) 

while  on  setting  k  =  0(tfg  =  <?z)  s 

q  =  (l  -  sin  <p)p  -  c  cos  <p  (19) 

Again  both  (18)  and  (19)  satisfy  the  hypothesis  of  Haar  and  v„  Karman,  intermedi¬ 
ate  formulae  being  obtained  by  using  other  values  of  k. 

The  ranges  of  possible  values  of  q  represented  by  (11)  and  (17)  are 
illustrated  in  Figure  5>  for  the  particular  case  c  *  1  and  q>  *  20°  . 

In  the  case  of  axial  compression,  the  usual  formula  (12),  represented 
by  the  upper  boundary  of  the  shaded  zone,  gives  the  highest  possible  values  of  q. 
Loading  paths  which  start  at  a  hydrostatic  pressure  state  (q  =  p)  must  cross  a 
zone  in  which,  according  to  limiting  equilibrium  theory,  failure  is  possible  with 
a  non-uniform  stress  distribution* 

In  the  above  illustration  of  alternative  stress  states,  it  is  assumed 
that  the  principal  stresses  are  in  constant  ratio,  although  there  is  no  obvious 
reason  why  this  should  be  so.  If  the  assumption  is  abandoned  and  k  made  a  func¬ 
tion  of  the  radius,  many  other  stress  distributions  all  satisfying  the  yield 
condition  and  the  boundary  conditions  become  possible. 

A  choice  between  all  these  solutions  can  be  made  by  introducing  a  stress- 
strain  relation  such  as  that  described  in  the  Introduction.  The  particular  stress- 
strain  relation  which  has  been  selected  narrows  the  choice  to  the  point  where  the 
external  tractions  become  uniquely  defined  and  a  rational  comparison  with  test  data 
becomes  possible.  With  the  aid  of  this  relation,  several  possible  velocity  fields 
will  be  found  which  are  associated  with  the  homogeneous  stress  state. 
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Plastic  Solution  for  Plane  Strain 


Suppose  a  single  plane  slip  band  of  width  t  forms  across  the  width  of 
the  specimen  (Fig.  6).  For  plane  strain  62  =  0  and  e^,  are  given  by 


i(ex  ♦  £y)  1  4  j  (6X  .  £y)2  +  Y. 
Substituting  in  either  (I4)  or  (5), 


'  \  2, 
2  I  2 


xy  t 

j 


e,r  =  y  tan  cp  .  (20) 

y  xy 


Thus  the  direction  of  motion  of  one  side  of  the  slip  band  relative  to  the  other 

subtends  an  angle  <p  with  the  direction  of  the  band.  The  corresponding  velocity 

field  is  ^ 

yu0  tan<P 


& 


yuo 

t 


(21) 


j 


in  the  slip  band  (t  y  ^0).  Both  in  the  case  of  axial  compression  (Fig.  6a) 
and  of  axial  extension  (Fig.  6b),  velocity  fields  have  now  been  found  that  are 
associated  with  the  (7)  and  (lU)  respectively.  The  remainder  of  the  specimen  is 
not  stressed  above  yield  and  the  solution  is  complete. 

Extent  of  the  deformable  region.  The  extent  of 'the  deformable  region 
can  now  be  found  by  using  the  theorem  due  to  Bishop,  Green  and  Hill ^ ^  which 
states  that  any  region  shown  to  be  necessarily  rigid  for  a  particular  stress  field 
(by  arguments  based  solely  on  the  geometrical  and  other  properties  of  that  field) 
must  be  rigid  in  all  complete  solutions.  The  slip  lines  form  the  characteristics 
both  for  the  stress  [12]  and  the  velocity  equations,  and  they  serve  to  isolate 
regions  at  the  ends  of  the  specimens  (as  indicated  in  Fig.  7)  which  are  neces¬ 
sarily  rigid  because  the  velocity  will  be  continuous  across  the  surface  of  the 
plattens.  Figure  7  can  represent  any  diametral  cross-section,  so  the  necessarily 
rigid  zones  are  conical. 
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Dilatation.  The  total  dilatation  can  be  computed  from  the  strains  (20) 
or  directly  from  energy  considerations.  By  adopting  the  latter  method,  the  dila¬ 
tation  rate  can  be  shown  to  be  independent  of  the  precise  nature  of  the  velocity 
field.  It  is  most  convenient  to  consider  separately  the  action  of  the  hydrostatic 
pressure  p  and  the  excess  axial  load  P.  Equating  the  work  done  by  the  external 
forces  during  unit  axial  compression  to  the  internal  energy  dissipation  computed 
from  (6), 

f 

P  -  p  A  =  c  cot  <P  *  6dv  =  c  A  cot  cp  (22) 

where  A  is  the  dilatation  of  the  entire  specimen  per  unit  axial  compression. 

In  the  case  of  failure  by  axial  compression,  (7)  applies,  and  it  can 


hence,  coirparing  (22)  and  (23), 


A  =  na2(N2-l) 


(2U) 


The  above  plane  strain  solution  establishes  a  possible  velocity  field 
for  the  homogeneous  stress  distribution  (obtained  by  setting  Oq  -  )  but, 

although  it  establishes  unique  external  pressure  values,  it  does  not  ensure 
uniqueness  of  the  stress  distribution,  due  to  the  presence  of  flats  on  the  yield 
surface. 
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Plastic  Solutions  for  Axial 


f!31 

The  analysis  given  below  follows  closely  that  given  by  Shield  J  for 
the  case  of  a  metal  ( <p  =0). 

Due  to  the  assumption  that  the  material  remains  isotropic,  the  principal 
directions  of  strain  will  coincide  with  the  principal  directions  of  stress]  hence 
Yrz  =  0,  or,  denoting  the  radial  and  axial  velocities  by  u  and  w  respectively. 


du  dw 

dz  dr 


*  0 


Other  strain  components  are 


3u  u  „  dw 

dr  *  e0  =  7  i  '  Tz 


(27) 


(28) 


Now  consider  possible  stress  states.  Fig.  U,  for  definiteness  setting  aa  “  ar  > 

cr,  *  oQ  and  o  **  cr  . 
b«  c  z 

Failure  with  axial  compression.  In  the  case  of  failure  by  axial  com¬ 
pression,  o^  >  cr^  and,  when  o0  is  the  intermediate  principal  stress,  only 
state  points  on  side  AF,  Fig.  U  need  be  considered.  On  side  AF,  except  at  A  and 
at  F,  6q  -  0  by  Table  I,  then  *  0  by  (28)  and  ez  =  0  by  (1*)  or  (5)j  hence 
no  solution  is  possible.  Possible  solutions  are  thus  immediately  restricted  to 
the  state  points  A  and  F. 

At  F  e q  is  positive  from  Table  I,  so  (it)  applies.  Recalling  (27)  and 
writing  (U)  in  terms  of  the  velocity  components  by  means  of  (28),  two  simultaneous 
partial  differential  equations  for  the  velocities  are  obtained: 


Shield  has  solved  these  equations  for  the  particular  case  N=1  and  solutions  for 
N  ^  1  can  be  found  along  exactly  the  same  lines.  A  first  deforming  zone  is  shown 
in  Figure  8a,  which  illustrates  one  half  of  the  cross-section.  Continuity  of 
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velocity  is  assumed  across  OB  so  that  u  =  w  =  0  on  OB  except  at  o.  Neither  (29) 
nor  the  boundary  conditions  involve  a  fundamental  length,  so  u  and  w  -will  be 
functions  of  t  only  (see  Fig.  8a),  and  (29)  reduces  to 


where  the  inverse  tangent  lies  between  0  and  n  .  The  strain  rates  are  in  the 
ratio 

er  :  6q  ez  =  tan2^  :  tan2(j|  +  -|)-tan2t  s  -1  (32) 

so  it  is  confirmed  that  er  and  6g  are  positive  when  is  negative.  The 
angle  t  varies  from  zero  to  (^  +  -|)  and  the  resulting  velocity  vectors  cover 
the  entire  range  of  permissible  directions  at  the  point  F,  as  indicated  in  Fig.  U» 
The  velocity  field  can  be  associated  with  one  and  only  one  stress  point,  so  that 
for  the  velocity  field  (31)  the  stresses  are  uniquely  defined  despite  the  presence 
of  flats  on  the  yield  surface.  As  a  consequence  it  is  not  necessary  to  investi¬ 
gate  state  point  A  as  a  possibility,  (if  this  is  done,  it  is  found  that  the 
strains  have  incorrect  signs.) 

The  total  dilatation  rate  has  been  found  previously  by  energy  consid¬ 
erations  j  it  can  also  be  found  by  integrating  the  sum  of  the  strain  rates  over 
the  deforming  region. 

The  nature  of  the  deformation  (31)  is  indicated  in  Figure  9a,  which 
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shows  the  deformed  shape  of  an  initially  square  grid  obtained  by  assuming  the 


initial  velocities  are  maintained  for  a  finite  deformation. 

Failure  with  axial  extension.  In  this  case  cl  <  cr_  and  failure  states  lie  on 
-  r  z 

side  CD  of  Figure  1*.  Except  at  C  and  D,  eQ  =  0  by  Table  I,  and  no  solution  is 
possible.  At  C,  eQ  is  negative,,  so  (5)  applies,  which  with  (2?)  gives 


3u  u  dw  - 

37  +  r  +  dz  tan2(£  >  f )  =  0 


3w 


(33) 


3u  3w 

'Si  +  37  “  0 

Adopting  the  deforming  region  indicated  in  Figure  8c  and  proceeding  exactly  as  in 
the  axial  compression  case,  the  velocities  are  found  to  be 

2 


u  =  -  -  7tan2({  -  f)  -  tan2  t 


(3U) 


w  =  -  —  tan"-*- y  tan2(^  -  •|)cot2\^  -  1 


where  again  the  inverse  tangent  lies  between  0  and  n  .  The  strain  rates  are  in 
the  ratio 

er  :  eQ  :  ez  =  -tan2t  s  -tan2(-^  -  4£)  +  tan2t  :  1 
which  is  consistent  only  with  the  state  point  C,  Figure  U.  Thus  the  homogeneous 
stress  state  is  the  only  possible  one  with  the  velocities  (3U).  The  nature  of  the 
deformation  is  indicated  in  Figure  9b. 

Other  solutions.  The  above  solution  serves  to  identify  the  homogeneous 
stress  state  with  an  axially  symmetric  mode  of  deformation.  That  the  mode  con¬ 
sidered  is  not  unique  is  established  by  the  alternative  solution  given  below. 

ri3i 

This  second  field  also  follows  closely  a  field  used  by  Shield1  J  for  the  case 
of  a  metal  (<p  »  0).  The  deforming  region  is  indicated  in  Figure  8b,  which  shows 
half  the  cross-section.  We  shall  consider  only  axial  compression,  for  which  the 
differential  equations  to  be  solved  are  (29).  In  the  region  OAB  the  field  used  is 


u 


2a 


tan  Or  +  |) 


w  =  tanOjj  -  f ) 


(35) 


1U 


which  satisfies  (29)  and  the  boundary  condition  w  -  1  on  AB.  The  velocity  across 
OA  is  proportional  to  the  distance  E  from  0*  so  it  is  assumed  that  in  AOC  the 


field  is  of  the  form 


u  ~  RF(t)  $  w  =  RG(t)  (36) 

where  F  and  G  are  functions  of  t  only.  The  condition  that  the  normal  velocity 
across  OA  be  continuous  (OA  is  a  slip  direction  for  the  stress  field  being  con¬ 
sidered)  requires 

F  tan(g  +  f )  +  G  =  [1  +  |tan2(|  +  |)  ]  /a  tan2(g  +  |)  +  1 

on  OA  and  F  »  G  =  0  on  OG,  Substitution  of  (35)  in  (29)  gives 
F(cost  +  sect)  -  F'sint  +  Gtan2(^  +  |)sint  | 


+  G'tan2(-J  +  |)cos  t  =  0  b 

i 

Fsint  +  F’cost  +  Gcost  -  G!sin  t  =  0  j 

✓ 

where  the  primes  indicate  differentiation  with  respect  to  t  •  The  variables  are 
separated  after  making  the  substitutions  F  =  A(f)eost  ;  G  =  B(t)sin  t  and  the 
solution  satisfying  the  boundary  conditions  on  F  and  G  is  found  to  be 


u=RF= 


w=RG= 


R  [tan2 (-^+^) cost  tan"^  y  tan2(-jj+^)cot2 1  -1  -  sint^tan2(j|+'|)-tan2 1] 

2 na  tan(g  +  |) 

Rtcost  V t an2  (j~+2  )  -tan2 1  -  sin  t  tan"1  \/tan2(^+2)cot^ t  -  1] 

na  tan(^  +  -|) 


where  the  inverse  tangent  lies  between  0  and  n <,  The  strain  rates  are  in  the 
ratio 


e  :  eQ  t  e  =  1  + 


1  -  T‘ 


cos“T  y 


-2  tan2(^  -  |) 


cos“-L  y 


where  y  =  tan(-£  +  -|)/tan  t  •  This  is  consistent  with  point  F  on  the  yield 
surface.  Figure  U® 
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All  the  solutions  given  above  are  associated  with  the  homogeneous  stress 
distribution,  the  axially  symmetric  modes  being  associated  with  that  distribution 
and  no  other.  While  it  cannot  be  proved  that  other  stress  distributions  are  im¬ 
possible,  in  view  of  the  flats  on  the  yield  surface,  they  appear  highly  unlikely. 
Moreover,  their  identification,  if  they  exist,  would  serve  no  useful  purpose 
because  the  present  solutions  give  unique  values  for  the  external  pressure  ratio 
and  the  dilatation  and  serve  to  show  that  both  plane  strain  and  axially  symmetric 
deformation  modes  do  in  fact  exist  for  an  ideally  plastic  material  that  fails 
according  to  the  Coulomb  yield  criterion. 

Experimental  Evidence 

The  above  analysis  enables  the  triaxial  test  to  be  used  to  check  the 
validity  of  the  hypothesis  of  Haar  and  v.  Karman  by  a  comparison  of  the  results 
of  extension  and  compression  tests.  The  complete  plastic  solutions  are  of  value 
in  that  greater  significance  is  now  placed  upon  the  Haar  and  v.  Karman  hypothesis, 
for  its  validity  becomes  a  necessary  (but  not  sufficient)  condition  for  the  soil 
to  behave  as  an  ideally  plastic  material. 

Several  series  of  tests  on  sands  under  extension  have  been  described  in 
the  literature,  and  the  results  are  summarized  in  Figure  10.  The  theoretical  lines 
in  this  Figure  were  computed  by  means  of  (l5)  and  (l6)  and  represent  the  extremes 
of  a/p  theoretically  possible  while  Og  remains  the  intermediate  principal  stress 
and  the  ratio  of  the  principal  stresses  is  constant.  The  test  results  are  plotted 
under  the  assumption  that  Oq  =  cr,  in  the  compression  test.  They  lie  in  the 
expected  zone,  but  there  is  a  considerable  discrepancy  between  the  results  ob¬ 
tained  by  Habib  3  in  1951  and  1953  and  those  obtained  by  Bishop  and 
Eldin^^  and  by  Kirkpatrick^  ^  in  1953  and  1957  respectively.  A  series  of 
tests  was  undertaken  in  an  attempt  to  obtain  more  information,  and  these  new 
results  are  also  indicated  in  Figure  10. 
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The  new  series  was  run  on  a  crushed  quartz  from  Lantern  Hill,  Mystic,  Conn., 
U.S.A.  The  quartz  was  almost  pure,  being  of  glass  making  quality.  The  fraction 
used  in  the  tests  was  that  retained  on  a  No.  $ 0  sieve  (opening  297  microns)  and 
passing  a  No.  U5  sieve  (opening  300  microns).  The  sand  was  placed  dry  by  tamping 
in  three  layers,  and  was  then  saturated  with  water  before  testing.  The  mean 
relative  density  achieved  by  this  process  was  70$. 

The  detailed  results  are  given  in  Table  II  and  in  Figure  11.  Two 
methods  of  reaching  the  failure  state  were  used  —  varying  the  axial  load  and 
reducing  the  lateral  pressure.  In  each  case  failure  was  deemed  to  have  occurred 
when  the  modulus  was  reduced  to  10  kg/cm^ ,  in  terms  of  whichever  pressure  was 
being  varied.  The  direction  of  approach  to  the  failure  point  is  indicated  in 
Figure  11  by  a  short  line  running  out  from  each  point.  A  mean  line  has  been 
drawn  through  the  compression  test  results,  and  the  angle  of  friction  computed 
from  this  line  (35.6°)  was  used  to  compute  the  two  theoretical  lines  shown  for 
the  extension  test  using  (18)  and  (19).  The  experimental  results  for  the  exten¬ 
sion  tests  fall  in  between  the  theoretical  lines,  the  tests  with  failure  by 
reducing  the  axial  pressure  showing  some  scatter  and  those  with  failure  by 
reducing  the  lateral  pressure  very  little*  Typical  points  are  also  plotted  in 
Figure  10  where  they  fall  close  to  those  established  by  Habib. 

Rates  of  volume  increase  per  unit  axial  compression  are  also  given  in 
Table  II.  These  are  the  average  values  at  failure.  It  should  be  noted  that  in 
series  C  and  D  the  specimens  were  extending,  so  negative  values  of  the  dilata¬ 
tion  per  unit  axial  compression  represent  an  expansion  of  the  specimen.  Theore¬ 
tical  values  of  dilatation  based  on  (2l|.)  for  compression  and  (26)  for  extension 
are  28  cc/cm.  and  -7.U  ec/cm.  respectively.  Thus  it  is  seen  that  the  actual 
dilatations  were  much  smaller  than  predicted  using  an  ideally  plastic  model.  No 
allowances  have  been  made  for  the  elastic  compression  of  the  specimens,  which 
would  probably  vary  according  to  the  intensity  of  the  hydrostatic  pressure,  and 
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would  in  all  cases  tend  to  decrease  the  dilatation * 

The  specimens  deformed  either  in  thin  slip  bands,  as  indicated  in  Figure  6, 
or  by  a  general  axially  symmetric  deformation,  as  in  Figure  9,  In  the  case  of 
failure  with  axial  compression,  it  was  usually  easy  to  decide  which  method  of 
failure  had  occurred  by  examining  the  specimen  at  the  conclusion  of  a  test.  In 
axial  extension  tests  it  was  much  more  difficult  to  decide,  a  typical  specimen 
being  that  shown  in  the  photograph  (Figure  12).  In  the  case  illustrated  an  axially 
symmetric  neck  has  formed  but  there  is  also  some  evidence  of  a  shear  plane  running 
from  top  left  to  bottom  right  of  the  sheared  zone.  In  all  the  extension  tests 
shear  was  confined  to  a  relatively  small  zone  similar  in  size  to  that  illustrated 
in  Figure  9b.  In  compression,  however,  shearing  frequently  extended  the  entire 
length  of  the  specimen  between  the  plattens.  This  result  is  as  expected:  changes 
in  geometry  tend  to  reduce  the  critical  cross-sections  in  extension,  leading  to 
increased  stresses  on  planes  which  have  slipped,  while  in  compression  the  cross- 
section  will  tend  to  be  increased  and  slip  may  well  be  transferred  to  a  smaller, 
and  hence  weaker,  plane. 

The  present  series  of  tests  tends  to  support  Habib's  conclusion  that  the 
effective  angle  of  friction  in  extension  tests  may  sometimes  differ  from  that  in 
compression.  Independent  tests  by  Bishop  and  Eldin  and  by  Kirkpatrick  show  that 
the  angle  is  identical  for  some  sands,  A  considerably  greater  volume  of  experimen¬ 
tal  evidence  is  of  course  needed  before  any  strong  assertions  can  be  made.  On 
the  evidence  to  date  it  does  appear,  however,  that  the  mechanical  behavior  of  sands 
may  vary  in  some  fundamental  manner,  governed  by  such  factors  as  grading,  particle 
shape,  porosity,  etc.  It  may  be  necessary  to  make  many  investigations  before  a 
clear  picture  emerges. 

Meanwhile  it  is  disturbing  to  find  some  sands  which  do  not  meet  the 
condition  of  yield  as  predicted  by  plastic  theory.  In  such  sands  the  standard 
limiting  equilbrium  approach  would  not  necessarily  result  in  conservative  esti¬ 
mates  of  carrying  capacity.  Previously,  there  was  a  distinct  possibility  that 
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all  limiting  equilibrium  states  would  prove  to  be  equal  to  or  to  underestimate  the 
actual  carrying  capacity*  as  if  the  soil  were  an  ideally  plastic  material. 

The  need  for  further  investigation  is  clear.  One  interpretation  of  the 
results  of  the  current  series  and  also  of  Habib’s  would  be  to  speculate  that  the 
intermediate  principal  stress  adopted  a  value  corresponding  to  k  =  0.4  in  (8). 

This  estimate  is  based  on  the  assumptions  that  the  principal  stresses  remain  in 
the  same  ratio  at  all  points  and  that  the  Coulomb  yield  criterion  applies.  Neither 
assumption  has  been  subjected  to  any  conclusive  empirical  test,  and  there  seems 
little  doubt  that  a  fuller  understanding  of  the  results  of  the  triaxial  test  must 
wait  upon  an  independent  investigation  of  the  influence  on  yield  of  the  intermedi¬ 
ate  principal  stress. 
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APPENDIX 


Uniqueness  of  applied  surface  forces  when  the  yield  surface  is  not  strictly  convex 
but  contains  flats. 

Theorem?-*  If  the  prescribed  surface  velocities  acting  on  a  body  composed  of  stable, 
ideally  plastic  material  are  uniform  in  magnitude  and  direction,  then  the  component 
of  the  resultant  surface  force  in  the  direction  of  motion  of  the  surface  is  uniquely 
defined. 


Proof  Consider  a  body  for  which  the  velocities  are  specified  over  the  region 
and  the  tractions  Ti  are  specified  over  the  region  of  the  surface.  Suppose 
a  solution  with  stresses  a.  strain  rates  e. .  and  resultant  surface  force  F-  on 
Au.  Let  etc.,  refer  to  any  equilibrium  state  that  does  not  violate  the 

yield  condition  and  satisfies  the  traction  condition  on  By  the  principle  of 

virtual  work 


r 


The  flow  rule  for  a  stable  material  (see  [93  )  requires  that 


This  theorem  has  been  proved  by  Hill  [l8]  for  the  more  restrictive  case  of  a 
strictly  convex  yield  surface. 
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Table  I 


a 


£b 


+  X2)N2 


-x: 


-  x: 


AB 


B 


X^N2 


-  X! 


x2n2 


■(xf  +  X2) 


x2n2 


BC 


-  x: 


x2n2 


-*s 


-  X! 


(X2  +  X2)N2 


CD 


x2n2 


.(x|  «■  X2) 


X2N2 


x2n2 


DE 


-  xu 


XpN2 


0 


E 


-  Xf 


2\m2 


(Xj  *  X|)H 


-  X, 


EF 


-  K 


X|N2 


X?N2 


-(x|  *  x|) 


1  FA 


X2N2 


~  Xz 


X 

1 


X  are  constants: 

6 
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Table  II 


I — 

Test 

Initial 

Porosity 

Lateral 
Pressure  p 

Mean  Axial 
Pressure  q 

Dilatation 

Rate 

% 

kg/ cm^ 

kg/ cm^ 

cc/cm 

A1 

U9 

0.6 

2.U 

U.9 

A2 

U9 

0.7 

2.7 

3.7 

A3 

51 

1.0 

3.7 

2.6 

Ah 

U7 

l.U 

5.U 

2.7 

A5 

5i 

1.6 

6.0 

2.0 

a6 

- 

2.1 

8.U 

a? 

U9 

2.6 

10.6 

2.2 

a8 

51 

2.6 

9.6 

1.2 

a  9 

U9 

2.6 

9.8 

1.6 

A10 

5o 

2.8 

10.7 

2.3 

All 

U9 

2.8 

10.5 

-» 

A12 

50 

3.0 

11.5 

1.8 

B1 

51 

0.6 

2.U 

5.3 

B2 

5o 

1.2 

U.U 

U.7 

B3 

U9 

1.8 

7.0 

«. 

Bh 

5o 

2.U 

8.U 

1.8 

Cl 

U9 

1.0 

0.3 

-l.U 

C2 

5o 

2.0 

0.6 

- 

C3 

5i 

3.0 

1.2 

0 

cu 

5o 

U.o 

1.3 

0 

c5 

U8 

5.o 

2.1 

0 

D1 

5i 

1.6 

0.6 

-0.3 

D2 

5o 

3.2 

1.2 

0 

D3 

5o 

3.7 

1.2 

-0.1 

DU 

5o 

U.7 

1.6 

0 

D  5 

5i 

5.6 

2.1 

— 

D6 

52 

6.U 

2.U 

+0.3 

Table  II  Results  of  triaxial  tests  on  crushed  quartz  (Lantern  Hill) 
Compression  tests 

Series  A  -  failure  by  increasing  axial  pressure 
Series  B  -  failure  by  decreasing  lateral  pressure 

Extension  tests 

Series  C  -  failure  by  decreasing  axial  pressure 
Series  D  -  failure  by  decreasing  lateral  pressure 


2h 


\  o 


ailures 


rigid  /  rigid 


FIG.8 


FIG.  10  Triaxial  extension  tests  on  sand.  Ratio  of  axial  to 
lateral  pressure  vs.  angle  of  friction  computed  from 
compression  tests. 


Kg /cm2 


FIG.  II  Triaxial  tests  on  crushed  quartz  sand 

(Lantern  Hill) 


FIG.  12  Appearance  of  a  specimen  after  an  extension 
failure,  showing  the  limited  extent  of  the  de* 
formed  zone.  (see:  FIG.  9b) 
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